PM/01-28 
TUM-HEP-434-01 
MPI-TH-32-01 
September 2001 



Loop induced Higgs and Z boson couplings to Neutralinos 
and implications for collider and Dark Matter searches 

A. Djouadi^ M. Drees^, P. Fileviez Perez^ and M. Muhlleitner^ 

^ Laboratoire de Physique Mathematique et Theorique, UMR5825-CNRS, 
Universite de Montpellier II, F-34095 Montpellier Cedex 5, France. 

"^Physik Department, Technische Universitdt Miinchen, 
James Franck Strasse, D-85748 Garching, Germany. 

^Max-Planck Institut fiir Physik (Werner-Heisenberg Institut) 
Fdhringer Ring 6, 80805 Miinchen, Germany. 



Abstract 

We calculate the one-loop induced couplings of two gaugino-like neutralinos to the Z 
and Higgs bosons in the Minimal Super symmetric Standard Model. These couplings, 
which vanish at the tree level, can be generated through loops involving fermions and 
sfermions. We show that, while the neutralino contribution to the invisible Z boson 
decay width remains small, the loop induced couplings to the lightest Higgs boson 
might be sufficiently large to yield a rate of invisible decays of this Higgs boson that 
should be detectable at future e^e~ colliders. We also study the implications of these 
couplings for direct searches of Dark Matter and show that they can modify appreciably 
the neutralino-nucleon elastic cross section for some parameter range. 
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1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) ^ is at the moment considered to 
be the most plausible extension of the Standard Model (SM). Out of the plethora of new 
particles contained in this model the lightest neutralino Xi plays a special role. It is often 
the lightest of all superparticles (LSP), which is absolutely stable if R-parity is conserved. 
This means that neutralino LSPs produced at colliders will be invisible, leading to the 
famous "missing [transverse] momentum" signatures for the production of superparticles. In 
particular, a final state consisting of two LSPs only would be entirely invisible, so that XiXi 
final states can contribute to the invisible width of the Z or neutral Higgs bosons. Moreover, 
if the LSP is stable, we expect some relic LSPs from the Big Bang era to still exist today, in 
addition to the well-known relic neutrinos and (microwave) photons. In fact, it was realized 
almost twenty years ago that x? is a good (cold) Dark Matter (DM) candidate 0. The 
LSP-nucleon scattering cross section, which determines the size of the expected signal in 
direct DM detection experiments 0, depends on the size of the LSP couplings to the Z and 
Higgs bosons; these couplings can also play a role in the calculation of the Xi I'slic density. A 
precise knowledge of these couplings is therefore important for both collider phenomenology 
and cosmology. 

The lightest neutralino x? couples to the Z boson and to the MSSM Higgs bosons only if 
it has a non-vanishing higgsino component. On the other hand, most models predict Xi to 
be dominantly a gaugino, in particular, a bino- or photino-like state. Moreover, an LSP with 
dominant higgsino component has a thermal relic density of the required magnitude only if 
its mass is in the TeV range, beyond the reach of near-future colliders and also beyond the 
range of masses that is usually considered to be natural. In contrast, a bino-like LSP is a 
good thermal DM candidate for masses in the 100 GeV range; masses in this range are more 
natural, and can be probed at present and near-future colliders. Since these gaugino-like 
states have suppressed tree-level couplings to gauge and Higgs bosons, loop contributions 
to these couphngs might be significant. 

In this paper we compute one-loop corrections to the neutralino couplings to the Z and 
Higgs bosons, where we only consider contributions with fermions and sfermions inside the 
loop [the other possible contributions, involving charginos, neutralinos and Higgs or gauge 
bosons, vanish in the bino-like limit for the x? ]• Our analytical results for these corrections, 
for arbitrary momenta and general fi — Jr mixing^, are given in Sec. 2. In Sec. 3 we apply 
these results to compute the loop-corrected invisible branching ratios of Z and Higgs bosons, 
and in Sec. 4 we study the loop-corrected LSP-nucleon scattering rate and relic density. 

For the neutralino couplings to Higgs and gauge bosons, we found in all cases corrections 
of up to a factor of two for reasonable values of the input parameters. The contribution 
of bino-like LSPs to the invisible decay width of the lightest MSSM Higgs boson might be 
measurable at future high-energy and high-luminosity e~^e~ colliders, while the contribution 
to the invisible decay width of the Z boson remains several orders of magnitude below the 

"'^Loop corrections to the Z exchange contribution to the bino annihilation cross section have previously 
been computed in Q. However, in this earlier paper /l — fn mixing has been ignored, and analytical results 
are only given for vanishing LSP three-momenta. Our numerical results agree qualitatively with their's. 
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present bound. On the other hand, the loop corrections can modify appreciably the LSP- 
nucleon scattering cross section in some MSSM parameter range. 



2 Neutralino couplings to Higgs and Z Bosons 
2.1 Tree- level couplings 

At the tree level, the couplings of the neutralinos to the Z boson are given by the vertex^ 



rJ(Zx°x;) = 7,75 [iV.3iV,3 - N,,N,,\ ^ 0,07,75, 



. 9 



2c 



(1) 



w 



while the neutralino couplings to the neutral CP-even Higgs bosons cf) = h,H and to the 
CP-odd boson A read |, 0]: 



(2) 



[{Ni2-tan9wNii){dANj3 + eANj4) + i ^ j]j5 = 9ax°x°^5- (3) 



Here, g = e/sw is the SU(2) gauge coupling with = 1 — c^r = ?>vo?9y^. The quantities 
ci<i,,e<i, ($ = h,H,A) are determined by the ratio tan/3 of the vacuum expectation values 
of the two doublet Higgs fields which are needed to break the electroweak symmetry in the 
MSSM, and the mixing angle a in the CP-even neutral Higgs sector: 



dn = —cos a , = sina , dA = sin/3, 
en = sina , Ch = cos a , ca = — cos/?. 

N is the matrix which diagonalizes the four dimensional neutralino mass matrix: 



(4) 



M 



N 



Ml 

Ma 

-mzswcp nizcwcp 

nizswsp -nizcwsp 



-rnzswcp mzSwSp 

mzcwcp -mzcwsp 
-/i 
-fi 



(5) 



where Mi and M2 are the SUSY breaking masses for the U(l)y and SU(2) ^ gauginos, is the 
higgsino mass parameter, and Sj3 = sin (3, etc. This matrix can be diagonalized analytically 
[^, but the expressions of the neutralino masses and the Nij matrix elements are rather 
involved. However, if the entries in the off-diagonal 2x2 submatrices in eq. are small 
compared to (differences of) the diagonal entries, one can expand the eigenvalues in powers 

^We assume that all soft breaking parameters as well as /i are real, i.e. conserve CP. We can then work 
with a real, orthogonal neutralino mixing matrix N if we allow the eigenvalues m-^o to be negative. Note 
also that the vertex factors F given in the text are 2i times the coefficients of the relevant terms in the 
interaction Lagrangian. 
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m^o ~ Ml —2 

mi 

rriyo ~ M2 -jK 



X4 



(Ml + /iS2/3) 



(M2 + /iS2/3) c^; 



"^1(1 - S2/3) 



'14/ 



;U + Mi /i + M2^ 
2 



/i-Mi /i-M2 



(6) 



In this analysis, we are interested in the situation > Mi, M2 and /i^ ^ m|. In this case 
the lighter two neutralinos will be gaugino-like. If |Mi| < IM2I, which is the case if gaugino 



masses unify at the same scale where the gauge couplings appear to meet the lightest 
state will be bino-like, and the next-to-lightest state will be wino-like. The two heaviest 
states will be dominated by their higgsino components.^ The eigenvectors of the mass matrix 
(PD can also be expanded in powers of mz- We find for the bino-like state [H, M: 



iVii 

iVl2 

iVii 

iVi3 

iVii 

iVi4 

iVii 



1 + (iVi2/iVii)' + (iVi3/iVii)2 + (iVi4/iVii)2' 



-1/2 



Ml Ml - M2 
■SffjU + c^Mi 

/i2 - Ml 



rrizsw 
-nizsw- 



+ 0{m%)- 



27) 



0{ml). 



(7) 



The corresponding expressions for the wino-like state read: 



A^22 

A^22 
A^23 
iV22 
iV24 

iV22 



1 + {N21/N, 



22 J 



(iV23/iV22)' + (iV24/iV22)' 



-1/2 



m?zSwCw S2I3IJ' + M2 



M| M2-M1 



Spfl + CpM2 2^ 

-mzcw — ^ TfT- + ^h^z) 



mzcw 



/i2 - Ml 
+ 5^M2 
/i2 - M| 



Note that the higgsino components of the gaugino-like states start at 0{mz), whereas the 
masses of these states deviate from their |yu| — > 00 limit (Mi and M2) only at (9(m|). 

Inserting eqs. (0) into eq. (|1|) one sees that the coupling of the LSP neutralinos to the Z 
boson only occurs at order m|, while eqs. (@) and (^ show that the LSP couplings to the 



^Loop corrections can significantly change the mass splitting between the two higgsino-like states |T3, |l3[. 
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Higgs bosons ^ — h, H, A already receive contributions at 0{m,z)- 



swmz 



(9) 



Similar expressions can be given for the couplings of the wino-like state. This suppression 
of the tree-level couphngs follows from the fact that, in the neutrahno sector, the Z only 
couples to two higgsino current cigenstates while a Higgs boson couples to one higgsino 

and one gaugino current eigenstate, together with the fact that mixing between current 
eigenstates is suppressed if ^ mz- Both kinds of couplings thus vanish as oo. 

The situation in the chargino sector is somewhat similar to what has been discussed so 
far. The diagonalization of the chargino mass matrix 



M 



c 



M2 \/2mwSp 



(10) 



leads, when expanded in powers of mw, to the two chargino masses: 



m^± -M2- (M2 + ixs2p) ~ m^o , m^± ~ /X + ^ + /x) , (11) 

so that for — > cxo, the lightest chargino corresponds to a pure wino state while the heavier 
chargino corresponds to a pure higgsino state. The couphngs of the neutral Higgs bosons to 
chargino pairs g'^{^xfj<^) are suppressed in this limit, and only the couplings g'°($xf X^) 
survive. Moreover, neither the W nor the charged Higgs bosons couple the LSP to the lighter 
chargino in this limit, g^{H^XiXi) ~ 0{mw/^J), g^iW^XiXi) ~ Oim^/ii^). 

However, the Z boson does have full-strength couplings to pairs of charginos xf X?i and 
LEP2 limits imply that even the lighter chargino is too heavy to be produced in the decay 
of the lighter neutral Higgs boson h. Moreover, the heavy Higgs bosons can always undergo 
unsuppressed decays into at least some SM fermions; in fact, decays involving 6— quarks are 
usually enhanced at large tan (3. We therefore do not expect heavy Higgs decays into modes 
that vanish for — > cxd to be significant even after loop corrections have been applied. For 
these reasons we will not discuss the chargino sector any further in this paper. 



2.2 One- loop induced couplings 

At the one-loop level, the couplings of the lightest neutralinos to the Z and Higgs bosons can 
be generated, in principle, by diagrams with the exchange of either sfermions and fermions, 
or of charginos or neutralinos together with gauge or Higgs bosons, in the loop. However, 
the latter class of diagrams can contribute to the couphngs of Higgs bosons to neutralinos 
only if one of the particles in the loop is a higgsino. Similarly, these gauge-Higgs loops will 
contribute to the coupling of bino-like LSPs to the Z boson only if at least one particle in 
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the loop is a higgsino-like state. These loop contributions will thus be suppressed by inverse 
powers of in addition to the usual loop suppression factor. We therefore do not expect 
them to be able to compete with the tree-level couplings that exist for finite \fi\, see eq. 

We thus only consider diagrams with fermions and sfermions in the loop, as shown 
in Fig. 1. In the case of the ZxiXi coupling, all three generations of matter particles 
will be involved since they have full gauge coupling strength. In the case of the 
couplings, only the third generation (s)particles, which have large Yukawa couplings, can 
give significant contributions to the amplitudes. Note that in the bino limit there is no wave 
function renormalization to perform, since the tree-level couplings are zero. Off-diagonal 
wave function renormalization diagrams could convert one of the gaugino-like neutralinos 
into a higgsino-like state, but this kind of contribution is again suppressed by and 
can thus not compete with the tree-level coupling. In case of the Z coupling, both external 
gauginos would have to be converted to higgsino-like states, which is obviously only possible 
at the two-loop level. 




Figure 1: The Feynman diagrams contributing to the one-loop couplings of the lightest 
neutralinos to the Z and ^ = h, H, A Higgs bosons. Diagrams with crossed neutralino lines 
have to be added. 



We have calculated the contributions of these diagrams for arbitrary momentum square 
of the Higgs and Z bosons, finite masses for the internal fermions and sfermions as well as for 
the external LSP neutralinos, and taking into account the full mixing in the sfermion sector. 
The amplitudes are ultra-violet finite as it should be. The contributions from diagrams a) 
and b) to the couplings are separately finite for each fermion species.0 The Zx^xl 

couplings are finite once contributions from both sfermion mass eigenstates of a given flavor 
and from both diagrams in Fig. 1 are added; diagrams a) and b) are separately finite after 

^Loops containing charged gauge bosons and gauginos only can contribute to the ZxiXi couphng if Xi 
is wino-like. However, these contributions are separately gauge invariant only if the Z boson is on-shell. 
Since in this scenario Tn^± ~ m-^o , LEP searches imply that xi is too heavy to be produced in the decays 
of on-shell Z bosons. If the Z boson is off-shell, box-diagrams with two charged gauge bosons in the loop 
have to be added to obtain a gauge-invariant result. 

^The contribution of diagram a) is finite only after summation over both sfermion mass eigenstates. 
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summation over a complete generation of (s)fermions. We have performed the calculation 
in the dimensional reduction scheme [jT4|; since the one-loop couplings are finite and do not 
require any renormalization, the same result is obtained using the dimensional regularization 
scheme once the contributions from a complete SU(2) doublet have been added. ^ The 
results are given below for a general gaugino ^ Mi,M2, for arbitrary ordering of Mi 
and M2). 

The one-loop induced Z boson vertex to a pair of lightest neutralinos is given by {pi is 
an incoming, p2 an outgoing momentum): 
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47r2 cw 



with the color factors A*"^^-^-*, and 

4/) = {af + VfC2ef)vi[-m}Co{h)-mloCo{f2)+^mloCt{f2) 

- 2(m|o +pi ■p2)C2+(/2) - 2(m|o -pi ■p2)C2-(/2) - 2C°(/2 
+ {aj - VfC2e^)vi [-mjCoih) - m|oCo(/i) + 4m|oC+(/i) 

- 2(m|o + pi ■ P2)C^2+(/i) - 2(m|o - ■ P2)C^2"(/i) - 2C2°(/i 
+ 2m;m^oa/S2e/t;i + v^) [Co{fi) - Coih) - 2C+(/i) + 2C+(/2) 



- {aj + Vf)v2m'j 
+ {af-Vf)v2 



C,{h)cl + C,{f2)sl 



m|oCo(/i) + 4m|„C+(/i) - 2C2°(/i) 



2(m|o +Pi ■P2)C^2^(/i) -2(m|o -pi •p2)C^2"(/i) 



+ st 



C20jVi 



V2 



, (-m|oCo(/2) + 4m|oC+(/2) - 2C2°(/2) 

2(m|o +Pl ■p2)C^2^(/2) -2(m|o -pi ■p2)C^2"(/2) 
fzfifi^2ifl^ fl) ~ /z/2/2^2(/2! /2) 
/z/i/iC'2 (/l) /l) + ^effzf2f2^2if2, f2 



+ S2e-{2vr + V2)fzjjC'2ifiJ~2) 



(12) 



(13) 



= (a/ + VfC2ej)'2m^oVi [Ctih) - 2C2 (h) 
+ (a/ - t;/C2ep2m^ot;i [C+{f\) - 2C,{f\) 
+ 2afS2ef{vi + t^a)^/ [-C+(/i) + C+(/2) 



®In the MS scheme the contribution from a given fermion species contains mass-independent terms 
proportional to some combinations of Z ff coupUngs. These constant terms add to zero when summed over 
a complete generation of fermions and sfcrmions. 
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-e-/z/i/i^2 /i) + ^e-fzf2f2^2 U2, 12) 



+ 2/z/i/2 [(^1 + V3)mfC^ (/i, /a) - 520^,(2^^1 + t^2)m^o<^2 /2) 

Here we have used the notation 



(14) 



V2 
V3 



{gQfNiitsiia9wf , vi 



2Qf 
-'3 



-Vo 



jf / 

-2 + -f + 2 1 
Qf \ 

Nl2 



'3 ) 



N- 



12 



N- 



12 



Qf J iVii tangly Qf vA^iitan6'vy 



21 



2Qf \Nuteiia9w 
The Z boson couphngs to sfermions are given by: 



(15) 



(16) 



The angle 6'^ is the mixing angle arising from the diagonalization of the sfermion mass 



matrices [compare eq. (0) below] and Sg. 



sin^ Or. Qf and // denote the electric 



charge and weak isospin of the fermion /, respectively. The Passarino-Veltman three-point 
functions iTBI, defined as 



^0,1,2 (/) = C'o;i,2 (? ,m^<i,mf,mj,mj); 
Co,i,f{fi, h) = C'o^2"(?^ "^|o, m|, m|, mj), 



(17) 



with q the momentum of the Higgs or Z boson, can be found in Ref. ||13|].fl 

In the same notation, the one-loop Higgs boson couplings to the LSP neutralinos in the 
bino limit are given by: 



47]-2 

9 
47r2 



f 



(/) 



75 



^Some care has to be taken if ^ or = 4m?o, since standard expressions for the foop functions con- 
tain spurious divergences in these kinematical situations, which are characteristic for LSP-nucleus scattering 
and LSP annihilation, respectively. This problem is discussed in the Appendix of Ref. p3| . 
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where 



{^2.,(^i + vs) [- {ml +m} + m|„) Co(/i) + 4m|oC+(A) 

+ {ml +m} + m|o) Coih) - 4m|oC+(/2)' 
+ 2(vi + V2cl_)mfm^o [Co(/i) - 2C+(/i) 

+ 2(i;i + V2sl_)mfm^o [Co(/2) - 2C+(/2) 
C^/i/i |-S2e^-K + V3)mfCo{fi, fi) + 2{vi + V2cl _)m^oCt {h, 
^0/2/2 {«2ej(^i + V3)mfCo{f2, h) + 2{vi + V2sl,)m^oC^{f2, h 
^4>hh {-2c2e^-(^;i + V3)mfCo{fi, f2) - 2s2e^-V2m^oC+{fi, /a)} ; 



(19) 



2mw 



Vl 



+ V3)s2e^- [{ml -m)- m|o) Co(/i) - (m| - - m|o) Co(/2) 



+ 2(111 + V2cl.)m^omfCo{fi) + 2{vi + V2SQ.)m^omfCo{f2)^ 
+ C^jj^ {-2mf{vi + V3)Co{fi, h) + 2m^o{2vi + V2)s2ef^{k h)} ■ 
The Higgs-fermion-fermion couphng constants are given by 

cos a — sin a 



sin [3 
sin a 



Qhuu 
9Huu n 

sm p 

QAuu = cot /3 



9hdd 
QHdd 



COS P ' 

COS a 



cos /3 ' 
gAdd = tan /3 , 



while the Higgs-sfermion-sfermion couphng constants read: 



ChU2U2 



mz 

Cw 

m,z 



Sf3+a 



TU 2 z 

-'3% - ^uSw(^2eii 



m^Qhuu muS2ea 



w 



1112 



c. 



hdidi 



c ~ ~ 

hd2 d2 



^ hd\ d2 



Cw 

mz 

Cw 

— -Sp+a Vt^l, ~ QdSwC2e^ 

Cw '- ''■ 

— S(3+a + QdSwC2e^ 

Cw '- 
mz 



mw 

i2. 



2m 



w 



UUl ) 



"'Z r™ 2 , ^ 2 1 f^ughuu . muS2eu r . , 1 
S(5+a \J-3 Sfil- + Vw'5w^C26(j — h [^uS'huu + l^gHuu] ] 



mw 



2m 



2m 



w 



[Aughuu + iJ-gH uu\ ^26ir '1 



mighdd ^4820^ 



[^dghdd — f^gHdd] 



Cw 



Sp+a [Qd^w ~ 2 



mw 
_ mg 

^ 2mw 



mw 2mw 

dghdd _|_ _ d \^^^g^^^ _ nQjj^^ ; 



2mvK 

[^dghdd — fJ'gHdd] C29/, 



m^gHuu 'muS2e^ 



(20) 



(21) 



mw 



2m 



w 



[^ugHuu ~ i^ghuu] ; 



C'HU2U2 



c 



Hdidi 



c 



Hd2d2 

^Hdid2 
^Adid2 



mz 

Cw 

mz 

Cw 

mz 

Cw 

mz 

Cw 

mz 



Cf3+a 



^3^ da + QuSwC29a 



m^QHuu m^S2B^ 

+ — V^uQHuu — f^dhv 



2m 



w 



^29,-,, 



mw 

— l^ugHuu — f^ghuu] C20fl ; 



2m 



w 



m^^gudd mdS2e^ 



mw 



2m 



w 



[^dgHdd + iJ'ghdd\ 



-^3 -^e, + QdSwC2e, 



m^gndd ^ _d_^ [Aagndd + fJ^ghdd] 



Qasl - 4/2 



Cw 

—^{Au cot/3 + /x); 
2mw 
md 

2mw 



mw 
md 



2m 



w 



2m 



w 



[AdgHdd + l^ghdd\ C26»j; 



(Arftan/? + /i). 



We use the following convention for the sfermion mass matrices: 



M) 



m'j + mf Af 

mj Af m) + 



m 



m 



with m^^ = m^^ + Q/s^ C2/3 



Af = A/ -/i(tan/5)-24 



(22) 



(23) 



They are diagonalized by 2 x 2 rotation matrices described by the angles 6f, which turn the 

current eigenstates, fi and /r, into the mass eigenstates /i and /2; the mixing angle and 
sfermion masses are then given by 



2mfAf 



2 2 I 

m? = mr H — 
/i,2 / ^ 2 



C2e, 



f m% — m% ' 



1 1 ' 



(24) 
(25) 



3 Higgs and Z Boson Decays 
3.1 Invisible decays of the Z boson 

The tree level and the one-loop induced axial vector couplings of the Z boson to pairs are 
shown in Fig. 2a for the input values tan /5 = 15 and /i = 1 TeV. We work in an unconstrained 
model with non-universal boundary conditions for the gaugino mass parameters at the high 
energy scale,^ and set the bino and wino mass parameters to M\ = 30 GeV and M2 = 120 
GeV, respectively. For definiteness we assume a common soft SUSY breaking scalar mass 
for the three generations of sleptons, m, 
squarks = m^ = mj 



El 

mg, with mq 



m 



mj, and a common mass term for the 
2mj. For the trilinear coupling, which will play 



®In models with universal gaugino masses at the GUT scale, the bound on the lightest chargino mass 
from negative searches at LEP2, m~± > 104 GeV |8|, will constrain a bino-like neutralino to have a 

mass larger than mz/2; the decay Z — > XiXi would then be kinematically closed. 
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a role mainly in the stop sector, we choose the value At = 2.9mg. The trilinear couplings in 
the (i— squark and slepton sectors, which are not important here, have been set to zero. 




Figure 2: The Z boson axial vector coupling (left) and branching ratio (right) to pairs of the 
lightest neutralino as functions of the common slepton mass. The dashed lines show the tree- 
level results; the solid line in the left frame shows the real part of the one-loop contribution 
eq. (|1^) to the axial vector coupling, while the solid line in the right frame shows the total 
loop-corrected branching ratio eq. (p6D. These results are given for tan/5 = 15,/i = 1 TeV, 
rrig = 2m;-, At = 2.9m^ and gaugino masses Mi = 30 GeV, M2 = 120 GeV. 



The tree-level coupling S'^^Oj^o is very small if Xi is bino-like, being of 0{m\/ ^^). As 
can be seen from Fig. 2a, it is less than 10~^ for the chosen parameters [and of course does 
not depend on slepton and squark masses]. The real part of the one-loop induced coupling, 
^^XiXi' ^^'^ on-shell Z bosons [i.e. with = m|] receives contributions from both massless 
and massive SM fermions and their superpartners. In the former case the real part has an 
"accidental" zero for sfermion masses slightly above the Z mass; this happens e.g. for the 
(s)neutrino contribution at m,> ~ 94 GeV.^ The real part of this contribution is positive 
(negative) for smaller (larger) sneutrino masses; its absolute value reaches a maximum at 
rrii) ~ 155 GeV. No such zero occurs for the imaginary parts. Moreover, the (s)top loop 
contribution is much larger than that from the other (s) quarks, even if the stop masses are 
not reduced significantly compared to the other squark masses. When combined with the 
"accidental" suppression of the (s)lepton contribution for slepton masses close to their present 
lower experimental bound, this implies that the dominant contribution to the real part of 
this coupling usually comes from (s)top loops. Of course, these loops cannot contribute to 

^This statement holds, strictly speaking, only in the MS scheme, where the contribution from each 
(s)fermion species decouples separately when — > 00. In the DR scheme decoupling only occurs after 
summation over a complete generation. 
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the imaginary part, which receives its dominant contribution from (s)lepton loops. For large 
masses, sparticles decouple and the loop correction to the coupling vanishes asymptotically. 

The partial decay width for the decay of a Z boson into a pair of LSP neutralinos is 
given by: 

r(Z^XiXi) = dzx^o+dzx^o , (26) 

where = (1 — 4m|o/m|)^/^ is the velocity of the neutralinos in the rest frame of the 
Z boson. In the usual perturbative expansion the one-loop correction to a decay width 
originates from the interference of tree-level diagrams with the real parts of the corresponding 
one-loop diagrams. However, in our case the tree-level decay width will vanish in the limit 
\n\ — > oo. In this limit the one-loop diagrams represent the leading order contributions. We 
therefore also include the squared one-loop contribution in eq. (^). However, for < 1 
TeV the numerically most important correction to the tree-level result usually still comes 
from the interference between the tree level and the real part of the one-loop couplings. 
Note finally that the derivative couphng g^^-o-o of eq. (IT^) does not contribute to the decay 
width of real Z bosons, since the product of the Z polarization vector with the sum of the 
outgoing LSP momenta vanishes. 

The Z decay branching ratio, BR{Z = r(Z ^ XiXDf^T with T*/ = 2.5 GeV, 

is shown in Fig. 2b. In the Born approximation, the branching ratio is very small, ~ 1 ■ 10"''. 
It can be enhanced by the loop contribution by a factor of > 2 in the low slepton mass range. 
However, the branching fraction is still well below the experimental limit on the invisible Z 



boson decay width |]17 



ATTUp ~ ±1.5 MeV, (27) 

which would require a branching ratio in the permille range for detection. One can find 
scenarios with somewhat larger loop corrections than shown in Fig. 2, if parameters are 
chosen such that all slepton masses are near 150 GeV. However, even in this case we are 
still more than three orders of magnitude below the experimental limit (p7|). Of course the 
Z partial decay width into Xi pairs can be enhanced to a detectable level for lower values of 
the parameter /i. However, in this case the one-loop corrections [which remain more or less 
the same] would be relatively less important than for large 



3.2 Invisible decays of the Higgs bosons 

The real part of the coupling of an on-shell lightest h boson to a Xi pair is displayed, at 
the Born and one-loop level, in the left-hand frames of Figs. 3. It is again shown as a 
function of the sfermion masses, for tan/? = 15, /i = 1 TeV and pseudoscalar mass input 
values Ma = 200 GeV (upper frames) and 1 TeV (lower frames). We see that the tree-level 
coupling g^j^o^^ is typically a factor of three larger than the corresponding coupling g^zx°x^ 
of the Z boson. This can be understood from eqs. (|1]), (0) and (j^). In the "decoupling" 
limit M\ ^ m| the neutral Higgs mixing angle a satisfies Sq, = — C/j, Cq, = sp. In this limit 
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Figure 3: The lightest h boson couphngs (left) and branching ratios (right) to pairs of the 
lightest neutralinos as functions of the common slepton mass. Most parameters are as in 
Fig. 2, and we took Ma = 200 GeV (top) and 1 TeV (bottom). 
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the ratio of tree-level ZxiXi ^XiXi couplings becomes — rriz cot 2(3/ {2 fi); for the given 
choice of parameters this amounts to about 0.3, in good agreement with Figs. 2 and 3. On 
the other hand, for small sfermion masses the one-loop corrections to the hx^Xi coupling 
exceed o^-o-o by more than an order of magnitude. The corrections to the Z coupling 
to bino-like neutralinos only involve electroweak gauge couplings. In contrast, top quarks 
couple with 0{1) Yukawa coupling to the Higgs bosons, and for the given choice of large 
l^fl the (dimensionful) htiti coupling significantly exceeds the ti mass. Moreover, due to 
— tji mixing the lighter t mass eigenstate is often not only lighter than the other squarks, 
but also lighter than the sleptons. As a result, the loop corrections are numerically even 
more important in case of the Higgs coupling, even though the tree-level contribution to 
this coupling is nominally only suppressed by one power of If 1 is large, as in the 
present example, the loop corrections to the hxiXi coupling can even exceed the tree-level 
contribution. 

The variation of the one-loop contribution to the coupling is again mostly due to the 
natural decrease with increasing masses of the sfermions running in the loop, which decouple 
when they are much heavier than the h boson. For rriq ~ 420 GeV [i.e. m^- ~ 210 GeV], m^^ is 
near its experimental lower bound of ~ 100 GeV, due to strong ti — Ir mixing. This implies 
that nii^ will grow faster than linearly with increasing niq, which explains the very rapid 
decrease of the loop corrections. However, there is also a variation of the tree-level coupling 
for Ma = 200 GeV which, at first sight, is astonishing. It is caused by the variation of the 
mixing angle a in the CP-even Higgs sector, and to a lesser extent by the variation of M/^, 
due to the strong dependence of crucial loop corrections in the CP-even Higgs sector on the 
stop masses.0 In fact, for the set of input parameters with Ma = 200 GeV at small slepton 
masses, we are in the regime where sin a, which appears in the hxiXi coupling [and which 
enters the hbb coupling as will be discussed later], varies very quickly. This "pathological" 
region, where the phenomenology of the MSSM Higgs bosons is drastically affected, has been 



discussed in several places in the literature [19 



The partial widths for the decays of CP-even Higgs bosons, (p = h, H , into pairs of 
identical neutralinos are given by 



~0~0^ W 



(28) 



The branching ratios for the decays of the lightest h boson are shown in the right-hand 
frames of Figs. 3 for the same choice of parameters previously discussed. They have been 
calculated by implementing the one-loop Higgs couplings to neutralinos in the Fortran code 



HDECAY [|2T| which calculates all possible decays of the MSSM Higgs bosons and where all 
important corrections in the Higgs sector, in both the spectrum and the various decay widths, 
are included. We see that the situation is completely different from the case of the decay 
Z — ^ XiXv the branching ratio BR(/i — »• XiXi) can already exceed the one permille level 
with tree-level couplings. This enhancement is not only due the larger Quxlxl couplings 

lOWe remind the reader that LEP Higgs searches would completely exclude the MSSM in the absence of 
these corrections. 
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Figure 4: The lightest h boson couphngs (left) and branching ratios (right) to pairs of the 
lightest neutralinos as functions of Mi. The parameters are as in Fig. 2 with mi = 250 GeV 
and Ma = 200 GeV. 

compared to g^^o^o, but also due to the much smaller h total decay width]^, r^°* ~ a few 
MeV, compared to F^* ~ 2.5 GeV. After including the one-loop corrections, the branching 
fraction BR(/i XiXi) can be enhanced to reach the level of a few percent. 

The branching ratio is especially enhanced if the usually dominant decay into bb pairs 
is suppressed, i.e. if | sina| is very small; recall that the hbb coupling is oc sin a/ cos/3. In 
our examples this happens for Ma = 200 GeV and m^- ~ 250 GeV. In this case g^xO^o is 
only about half as large as in the decoupling limit Ma —>■ oo, but the loop contribution to 
this coupling is still sizable for this value of the sfermion masses, and has the same sign 
as the tree-level coupling, leading to a quite large total coupling. The branching ratio 
falls off quickly for smaller sfermion masses, since here sin a becomes sizable (and positive). 
Moreover, for ~ 210 GeV the tree-level and one-loop contributions to the couplings have 
opposite sign. The branching ratio also decreases when m^- is raised above 250 GeV, albeit 
somewhat more slowly; here the rapid decrease of the loop contribution is compensated by 
the increase of the tree-level coupling, which however does not suffice to compensate the 
simultaneous increase of sin^ a. 

Fig. 4 shows the dependence of the hxiXi coupling and of the corresponding h branching 
ratio on the mass of the LSP, m^o ~ Mi. We see that the tree-level contribution to this 
coupling depends essentially linearly on the LSP mass. Eqs. (g) and (|^) show that, for the 
given scenario where ~ ~ 1, this linear dependence on Mi originates from where 

^^For the choice of parameters in Fig. 3 the total decay width of h is similar to, or even smaller than, that 
of the Standard Model Higgs boson with equal mass. 
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the contribution with Mi in the numerator is enhanced by a factor of tan/3 relative to the 
contribution with /i in the numerator. Therefore the contribution oc Mi is not neghgible 
even though in Fig. 4 we have Mi -C On the other hand, the one-loop contribution 
to this coupling depends only very weakly on Mi. The small increase of this contribution 
shown in Fig. 4 is mostly due to the explicit m^o dependence of the loop coupling (|T9|); the 
change of A^i2 with increasing Mi, as described by eq. plays a less important role. The 
increase of the total coupling with increasing Mi nevertheless remains significant. However, 
the right panel in Fig. 4 shows that for m^o > 15 GeV this increase of the coupling is over- 
compensated by the decrease of the f3^ threshold factor in the expression ( pSj) for the h 
partial width. 

Once one-loop corrections are included, for certain values of the MSSM parameters the 
branching ratio for invisible h boson decays can thus reach the level of several percent even 
if Xi is an almost purely bino. This would make the detection of these decays possible at 
the next generation of e~^e~ linear colliders. At such a collider it will be possible to isolate 
e~^e~ — > Zh production followed hj Z decays {£ = e or /i) independent of the h decay 

mode, simply by studying the distribution of the mass recoiling against the pair. This 
allows accurate measurements of the various h decay branching ratios, including the one for 
invisible decays, with an error that is essentially determined by the available statistics P2|. 



Since a collider operating at ^/s ~ 300 to 500 GeV should produce ~ 10^ Zh pairs per year 
if I sin(Q; — /3)| ^ 1 one should be able to measure an invisible branching ratio of about 3% 
with a relative statistical uncertainty of about 2%. 

We now turn to the heavier MSSM Higgs bosons H and A. The couplings to the lightest 
neutralinos are shown in Fig. 5 for the same input parameters as in Fig. 2. As can be seen, up 
to a relative minus sign, the tree-level couplings of these two Higgs bosons are approximately 
the same since we are in the decoupling regime where c?^ ~ —dn, and ca — en with le^l <^ 1, 
see eq. (H). Eqs. (IH^) and (0) also show that the tree-level couplings of the heavy Higgs 
bosons exceed that of the light Higgs boson /i by a factor tan/?/2 [ignoring contributions 
to eqs. with Mi in the numerator]. On the other hand, the loop corrections are smaller 
in case of the heavy Higgs bosons. The corrections to the i^XiXi coupling are reduced by 
about a factor of 2 compared to the corrections to the /^XiXi coupling, mostly due to the 
relatively smaller coupling to ti pairs, see eqs. (^21) . The corrections to the Ax^Xi coupling 
are even smaller, since the CP-odd Higgs boson A cannot couple to two identical squarks. 
The contribution with two ti squarks and one top quark in the loop, which dominates the 
corrections to the couplings of the CP-even Higgs bosons for small m^-, does therefore not 
exist in case of the A boson. As a result, the corrections to the coupling of A are not only 
smaller, but also depend less strongly on m^; recall that for our choice of parameters m^^ 
increases very quickly as m^- is increased from its lowest allowed value of ~ 210 GeV, which 
comes from the requirement m^^ > 100 GeV. Note also that the Htt and Att couplings 
are suppressed by a factor cot (3 relative to the htt couphng, see eqs. ( pi]) ; this becomes 
important for large squark masses, where the contributions from Fig. lb are relatively less 
important. Altogether we thus see that the one-loop corrections are much less important for 
the heavy Higgs bosons. Note also that for the given set of parameters they tend to reduce 
the absolute size of these couplings. 
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Figure 5: The heavier CP-even Higgs boson H (upper plots) and the pseudoscalar A boson 
(lower plots) couplings (left) and branching ratios (right) to pairs of the hghtest neutralinos 
as functions of the common slepton mass. The parameters are as in Fig. 2. 
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The partial decay width r{H —>■ XiXi) of ^^e CP-even Higgs boson H is given by 
eq. (|28D with (p = H. Due to the different CP nature of the pseudoscalar Higgs boson A, the 
expression for its partial decay width differs slightly; it is given by: 

r(A^iK) = ^|,»-.,, + ,i-.,,f. (29) 

Again, because in the decoupling regime the CP-even H boson and the pseudoscalar A boson 
have almost the same couplings to Standard Model particles and to the neutralinos [at the 
tree level], their branching ratios are approximately the same. The one-loop contributions 
decrease the branching ratios by at most ~ 10 to 40%. Note that the total decay widths 
of the A and H bosons are strongly enhanced by tan^/5 factors [T{H, A — > 66) oc mlg\j^fj,j]. 
This over-compensates the increase of their couplings to neutralinos, so that their branching 
ratios into Xi pairs are far smaller than that of the light Higgs boson h, remaining below 
the 1 permille level over the entire parameter range shown. Moreover, the cross section 
for the production of heavy Higgs bosons at e^e~ colliders is dominated by associated HA 
production, which has a much less clean signature than Zh production does. Branching 
ratios of the size shown in Fig. 5 will therefore not be measurable at e^e~ colliders. In fact, 
they will probably even be difficult to measure at a fi^fi~ collider "Higgs factory"; recall 
that the Z factories LEP and SLC "only" determined the invisible decay width of the Z 
boson to ~ 0.1%, see eq. (p^). 



4 SUSY Dark Matter 

It is well known that in the MSSM with exact R-parity, the lightest neutralino is a good 
cold Dark Matter (DM) candidate g i. For a very reasonable range of supersymmetric 
parameters the relic density of a bino-like neutralino satisfies 0.1 < Q^oh"^ < 0.3, which is 



the currently preferred range p3[; here fl denotes the mass density in units of the critical 
density, and h is today's Hubble constant in units of 100 km/(sec-Mpc). Note that this 
argument singles out a bino-like LSP. A higgsino- or wino-like LSP would have a thermal 
relic density in this range only if its mass is around 1 TeV; such a large mass for the lightest 
superparticle would be difficult to reconcile with finetuning or naturalness arguments. 

The two relic neutralino search strategies that suffer least from our lack of knowledge of 
the DM distribution throughout our galaxy are "direct" detection, where one looks for the 
elastic scattering of ambient neutralinos off nuclei in a detector; and the "indirect" search 
for high-energy neutrinos originating from XiXi annihilation in the center of the Earth or 
Sun. In both cases the signal rate is directly proportional to the LSP-nucleon cross section 
. This cross section in turn depends on the neutralino-quark interaction strength, on the 
distribution of quarks in the nucleon, and on the distribution of nucleons in the nucleus. 
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Figure 6: Diagrams contributing to the effective neutralino-quark interactions. 



Three classes of diagrams contribute to the neutrahno quark interaction: the exchange of 
a Z or Higgs boson in the t— channel, and squark exchange in the s— or u— channel, see Fig. 6. 
Z exchange only leads to a spin-dependent interaction, while Higgs exchange contributes 
only to the spin-independent interaction, and squark exchange gives both spin-dependent 
and spin-independent contributions, where the latter arises only due to combinations of 
couplings that violate chirality; see Ref. for details. In general the spin-independent 
contribution is more important since it leads to a coherent coupling to heavy nuclei]^ and 
the Higgs exchange contribution is usually bigger than the squark exchange contribution. 

We saw in Sec. 2.1 that in the pure gaugino limit the XiXiHiggs couplings vanish at the 
tree level. We just stated that the most important contribution to the spin-independent 
neutralino quark interaction often comes from Higgs exchange. In this Section we therefore 
study the effect of the fermion-sfermion loop corrections to the Higgs couplings to neu- 
tralinos on the neutralino-proton cross section. We use the expressions in Ref. these 
expressions also include contributions from effective LSP-gluon interactions. We also take 
into account leading SUSY QCD corrections to the scattering cross section ||25|, and assume 
the "standard" value for the strange contribution to the nucleon mass, ms(p|ss|p) = 130 
MeV. 

As in Sec. 2 we use non-universal soft breaking terms but keep all squark masses identical. 
All experimental constraints on sparticle and Higgs boson masses are taken into account. 
We show two illustrative examples. In Fig. 7 we take m^o ~ M\ = 50 GeV, and vary fi from 
700 to 1050 GeV. The values of the other parameters are similar to those of Fig. 2. The 
upper panel shows the predicted relic LSP proton scattering cross section, while the lower 
panel shows the thermal LSP relic density Vt^oh?. In the case at hand this latter quantity 
is essentially determined |p by Xi pair annihilation into charged lepton pairs. The rapid 
decrease with increasing /i is due to the decrease of m^^ , which is caused by enhanced — tr 
mixing; this not only trivially increases the fi propagator, but also allows Xi annihilation 
to proceed from an S-wave initial state 0. We include this figure here in order to illustrate 
that our choice of parameters leads to an LSP relic density of the required magnitude. 
Of course, this is not difficult to arrange, since the loop corrections we are interested in 

^^The x\ capture rate in the Sun, which contains very few heavy nuclei, can also get important contribu- 
tions from spin-dependent interactions. 
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Figure 7: The predicted neutralino-proton scattering cross section and the thermal rehc 

density (bottom), when one-loop corrections to the couplings of the Higgs boson to the 
LSP arc included (solid curves) or omitted (dashed curves); in the lower frame the two 
curves practically lie on top of each other. The values of the relevant parameters are: 
m^o ~ Ml = 50 GeV, = 400 GeV, Ma = 200 GeV, M2 = 105 GeV, = 2.73mg-, tan /3 = 
21,'me« = 210 GeV, and mo = 200 GeV. 
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are almost independent of soft breaking parameters in the slepton sector, which therefore 
always can be chosen to produce the desired relic density, as long as the LSP is not too heavy. 
Not surprisingly, the prediction of ^^oh? is practically not affected by the loop corrections 
presented in Sec. 2. 

On the other hand the effect of these corrections on the LSP-nucleon scattering cross 
section can be significant, as shown in the upper frame of Fig. 7, where they reduce the 
predicted scattering rate by up to a factor of 1.5. For the given choice of parameters the 
tree-level and one-loop contributions to the hxiXi coupling have opposite sign, as in the 
upper frames of Fig. 3; the absolute value of the sum of these contributions only amounts 
to about 20% of the absolute value of the tree-level coupling. However, the scattering cross 
section is actually dominated by the exchange of the heavier H boson. For the given choice 
of parameters its mass exceeds only by a factor of ~ 1.5, and its couplings to down- 
type quarks are enhanced by a factor ~ tan/5; moreover, we saw in Sec. 3 that 

exceeds (7?,-o,-o by a factor ~ tan/9/2. The absolute value of the ifY?y? coupling is also 
reduced here, as in Fig. 5, but only by ~ 10%. Another significant contribution comes 
from s squark exchange, since sl — sr mixing is enhanced for large ■ tan/3; this mixing 
increases the spin-independent part of the squark exchange contribution Of course 

this latter contribution, which interferes constructively with the if— exchange contribution, 
is not affected by the loop corrections of Sec. 2. Hence in the given case the net effect of 
these corrections is significantly smaller than in case of the invisible width of the h boson. 
Finally, the overall decline of the predicted LSP scattering rate with increasing fi is due to 
the reduction of the higgsino components of the LSP, see eqs. (|^), which leads to smaller 
tree-level LSP couplings to the Higgs bosons. Since the one-loop corrections to the Higgs 
couplings depend only rather weakly on n they become relatively more important when 
is increased. 

As a second example, shown in Fig. 8, we take m^o ~ Mi = 100 GeV, and slightly 
reduced values for and tan p. The reduction of tan $ implies less f^r/j mixing, so that 
the relic density is much less sensitive to n than in the previous example. On the other hand, 
the increase of the LSP mass while keeping the slepton masses essentially the same leads 
to an increase of the annihilation cross section, so that we now obtain a relic density of the 
required size in the entire range of n shown in this figure. The effect of the loop corrections 
to the Higgs XiXi vertices is somewhat larger here than in the previous example. Due to the 
reduced value of Ma we are now in a regime with strong Higgs mixing, i.e. | sina| is much 
larger than the value cos/3 it takes in the limit Ma — >■ oo. As a result the exchange of the 
lighter h boson is now more important than before. The loop corrections reduce the absolute 
value of the /^XiXi coupling by about a factor of two, which leads to a similar reduction of 
the predicted LSP scattering rate. 

We already saw in Sec. 3.2 that loop corrections to the HiggsXiXi couplings are maximal 
if ti is light and the Higgstiti couplings are large. This is again illustrated in Fig. 9, where 
we show the dependence of a ~o on At (top) and tan/3 (bottom). Reducing At from its 

upper bound (set by the experimental lower bound on uih) increases m^^ and reduces the 
hiiii coupling, which leads to a rapid decrease of the loop correction to cr~o_ . The At 
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dependence of the tree-level prediction is due to the t — i loop corrections to the MSSM 
Higgs sector, which are always included. This dependence is quite mild, except near the 
upper bound on At where both rrih and the mixing angle a depend sensitively on A^. On the 
other hand, the tree-level prediction for the scattering cross section depends quite strongly on 
tan P, due to the 1/ cos /5— behavior of the Yukawa couplings of down-type quarks. However, 
the negative loop corrections also increase in size with increasing tan/5, largely due to the 
tan/3 dependence of the Higgs couplings to t and ii (recall that we are far away from 
the "decoupling limit"). We thus see that in the case at hand the parameter dependence 
actually becomes weaker once loop corrections are included; however, this is not always the 
case. Finally, we note that scenarios with m^^ ~ m/j, where our loop corrections will be 
sizable, can also be found in models with universal soft breaking terms at the scale of 
Grand Unification, if the trilinear soft breaking parameter is sizable already at this scale; 
the condition ^ Mi is almost always satisfied in these models. 

In these two examples the corrections to the Higgs couplings to the LSP reduce the 
predicted scattering rate. It is easy to construct scenarios with even larger but positive 
corrections, by taking very large values for first and second generation squark masses as 
well as for Ma, while keeping third generation squark masses relatively small; this kind 
of spectrum is e.g. expected in "more minimal" supersymmetric models [^|. The largest 



contribution to the LSP-nucleon scattering cross section then comes from h exchange, so 
the corrections to this cross section would scale essentially like the square of the correction 
to the hxixl coupling discussed in Sec. 2.2. However, the predicted counting rate in such a 
scenario would be very low, well below the sensitivity of near future experiments. We finally 
note that the corrections to the ZxiXi couplings will have some effect on the spin-dependent 
LSP-nucleon scattering cross section. However, this cross section is usually dominated by 
first and second generation squark exchange, which in this case does not require violation of 
chirality. We therefore expect the loop corrections to the Z coupling to be significant only 
in scenarios where | Mi | ^ | /i | < m - j - but some sfermion masses are significantly smaller 
than 

It should be noted that currently the XiP scattering cross section can only be predicted 
to within a factor of 2 or so even if the fundamental LSP-quark interactions were known 
exactly, the main uncertainty coming from the poorly known size of the strange matrix 
element ms{p\ss\p). For the cases considered in this paper, the loop corrections calculated 
here can therefore shift the cross section by at most one "theoretical standard deviation". 
However, this theoretical uncertainty will presumably be reduced in future. We believe that 
the calculation presented here, together with the results of refs.|l^ and [^, reduces the 



theoretical uncertainty of the prediction of the hard scattering cross section for given SUSY 
parameters to the level of 10%. 



5 Summary and conclusions 

In this paper we have studied quantum corrections to the couplings of the Z and neutral Higgs 
bosons to gaugino-like neutralinos in the MSSM. We focused on the phenomenologically 
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most interesting case of a bino-like lightest neutralino x° as LSP, but our analytical results 
of Sec. 2 are valid for a more general gaugino-like neutralino, irrespective of whether it is the 
LSP. We found that these corrections can completely dominate the tree-level contribution to 
the coupling of the lightest CP-even Higgs boson. The corrections to the couplings of the Z 
and heavy CP-even Higgs boson are somewhat less significant, but can still amount to about 
a factor of 2. Since the CP-odd Higgs boson cannot couple to two identical sfermions the 
corrections are suppressed in this case. In all cases the corrections can be significant only if 
some sfermion masses are considerably smaller than the supersymmetric higgsino mass 
Both (s)lepton and (s)top loop contributions to the Z coupling can be important, since the 
experimental bounds for mi^ and are still quite close to Mz- The Higgs couplings receive 
their potentially largest corrections from loops involving third generation quarks and their 
superpartners; in the latter case the corrections are also quite sensitive to the size of the 
trilinear soft breaking parameter At [and Af,, if tan/? 3> 1]. 

Turning to applications of these calculations, we found that these corrections might 
change the predicted detection rate of Dark Matter LSPs by up to a factor of two even for 
scenarios where the rate is close to the sensitivity of the next round of direct Dark Matter 
detection experiments. The possible impact on the invisible width of the hghtest CP-even 
Higgs boson h is even more dramatic: it could be enhanced to a level that should be easily 
measurable at future high-energy e^e~ colliders, even if the LSP is an almost perfect bino. 
This would open a new window for testing the MSSM at the quantum level. 
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